Impact of climate change on river water temperature and dissolved oxygen: Indian riverine thermal regimes

The impact of climate change on the oxygen saturation content of the world’s surface waters is a significant topic for future water quality in a warming environment. While increasing river water temperatures (RWTs) with climate change signals have been the subject of several recent research, how climate change affects Dissolved Oxygen (DO) saturation levels have not been intensively studied. This study examined the direct effect of rising RWTs on saturated DO concentrations. For this, a hybrid deep learning model using Long Short-Term Memory integrated with k-nearest neighbor bootstrap resampling algorithm is developed for RWT prediction addressing sparse spatiotemporal RWT data for seven major polluted river catchments of India at a monthly scale. The summer RWT increase for Tunga-Bhadra, Sabarmati, Musi, Ganga, and Narmada basins are predicted as 3.1, 3.8, 5.8, 7.3, 7.8 °C, respectively, for 2071–2100 with ensemble of NASA Earth Exchange Global Daily Downscaled Projections of air temperature with Representative Concentration Pathway 8.5 scenario. The RWT increases up to7 °C for summer, reaching close to 35 °C, and decreases DO saturation capacity by 2–12% for 2071–2100. Overall, for every 1 °C RWT increase, there will be about 2.3% decrease in DO saturation level concentrations over Indian catchments under climate signals.

River water quality parameters such as River Water Temperature (RWT), and Dissolved Oxygen (DO) forms vital signs for defining the health of a river water body's ecosystem 1 . Global warming climates have also shown an adverse impact on RWT under intensification various climatological defining variables, majorly Air Temperature (AT) [2][3][4] . Intensification of RWT will have adverse impacts in terms of a decrease of river DO saturation levels, where most of the river water quality standards are defined based on such saturation levels 3 . Precisely, saturation DO is a prominent indicator of river water quality and is considered a standard measure to define the pollutant extent 5 . The influence of climate change on DO in relation to RWT can lead to water quality degradation and ecological distortion [6][7][8][9][10][11] . RWT is inversely related to DO concentration that every change in RWT affects the river's ability to self-purify by lowering the amount of oxygen that can be dissolved and utilized for biodegradation [12][13][14] . Hence, climate change impacts on RWT and saturation oxygen content are prominent in understanding the projected river water quality and possible alterations in quality standards under climate change warming signals.
Water quality modeling studies predicted depletion of DO under streamflow, RWT, and land use changes for various basins globally 1,7,[15][16][17][18][19] . Such studies modeled RWT based on regression models 18 and process-based stream temperature models 15,16 and river water quality models such as QUAL2K 17,19 . However, such studies are basin or river stretch specific, data intensive, and limits application for data sparse and ungauged locations with an emphasis on simulated DO levels in response to streamflow, RWT, and land use [16][17][18][19] . However, DO saturation level, which serves as a baseline to measure oxygen-based water quality by determining the oxygen concentration of unpolluted water depending on RWT, salinity, and oxygen partial pressure 1 and is prominent in defining the maximum permissible limits and standards for various river usages 20,21 , has not been assessed under climate change. Specifically, while some recent studies have looked at how climate change affects RWTs, the question of how climate change affects saturation DO have yet to be answered. More specifically, the direct integration of RWT predictions in the assessment of DO saturation concentration levels under climate change signals has not been quantified. Therefore, the present study aims to quantify the projected changes in DO saturation levels under RWT projections using the state-of-the-art Global Climate Model (GCM) projections. Furthermore, saturation DO is generally considered a desirable level of DO by the Pollution Control Boards (PCBs) in Waste www.nature.com/scientificreports/ Load Allocation Models (WLAM) for river water quality management 22 . Therefore, the study of climate change impacts on saturation DO levels can provide prominent insights for defining/alterations of the water quality standards under climate signals. Climate change has been demonstrated to have an impact on the relationship between RWT and DO concentrations in tropical rivers 7 . Tropical rivers receive more solar radiation and have higher RWTs 23 . For example, Indian tropical river systems experience the highest RWTs during low flow periods of non-monsoon and summer months 19,24 . Seasonality plays a vital role in the Indian river systems as maintaining flows in the summer season is a challenge leading to water quality deterioration. To this end, the assessment of DO saturation rates with respect to RWT is of much relevance for Indian river systems due to minimum flows and higher temperatures during non-monsoon seasons.
Accurate estimation of RWT is prominent and can be estimated based on thermal advection-dispersion models 25 , equilibrium temperature-based models 26 , statistical or machine learning (ML) models 27 , and hybrid models 28 . Unlike process-based models, ML models do not require many input variables, which are unavailable for many ungauged river systems and have been widely used as robust in RWT modeling in recent years 29 . In this context, regression models 18,27,[30][31][32][33][34] , classical ML models 35,36 , Artificial Neural Networks (ANN) [37][38][39][40][41][42][43][44][45] , has proven to be a viable technique for RWT forecasting. Most Indian River systems are burdened with data limitations and form a significant challenge for implementing process-based RWT models and promotes to implement regression or ML based approaches to predict RWT 35 . Therefore, given the limitations over data availability for Indian river systems, the present study stresses on use of ML based prediction algorithms that can address the data availability limitations. In this study, we used Long short-term memory (LSTM) model coupling with the k-nearest neighbor (k-NN) bootstrap resampling simulation technique (kNN-LSTM) to achieve a better prediction of RWT under data limitations for seven major Indian catchments with monthly RWT data. In summary, the objective of the study is to calculate the impacts of climate change on riverine thermal processes in India and possible variability in DO saturation levels with respect to RWT by using the kNN-LSTM model addressing sparse spatiotemporal RWT data forced with state-of-the-art climate change projections. The study evaluated the effect of climate change on DO saturation with respect to seasonal RWTs with an ensemble of 21 General Circulation Models (GCMs) using Representative Concentration Pathway (RCP) 8.5 scenario dataset output downscaled from the National Aeronautics Space Administration (NASA) Earth Exchange Global Daily Downscaled Projections (NEX-GDDP) dataset. The present study considered seven majorly polluted catchments of India 46,47 with various physiographic features to analyze climate change impacts on saturated DO with respect to predicted RWT using kNN-LSTM based ML model using NEX-GDDP projections.
Methods kNN-LSTM model. The present study considered the most widely known RNN architecture of LSTM to predict the RWT due to the superiority of using backpropagation through time and overcoming the vanishing gradient problem, and capable of learning long-term dependencies [48][49][50][51] . The LSTM consists of different memory blocks called cells. Each memory cell has an input gate, an output gate, and an internal state that feeds back into itself unaffected over time steps, which learns when it's time to forget about prior hidden states when to update hidden states given new data and be used to learn complex temporal sequences. These memories in LSTMs are called cells. This study used the ensemble of k-NN bootstrap resampling algorithm to simulate the data from historical records based on Raseman et al. 52 and LSTM model (kNN-LSTM) for monthly RWT prediction at the seven catchment sites of India with sufficient tests of performance measures of a model. For future RWT projections, RCP 8.5 scenario down-scaled projections of AT data were fed into the kNN-LSTM monthly prediction model. To train the kNN-LSTM model, the current AT and previous month time-lag of both AT and water temperature of the k-NN algorithm-based data as predictors for seven catchments of India at monthly timescale. The first month's water temperature is calculated based on the catchment mean from the historical record. The prediction of subsequent months proceeds as follows: where T w t+1 is the future RWT prediction at time t + 1 month; f is a non-linear function which is generated by the kNN-LSTM monthly model; T a t+1 is the future AT at time t + 1 month; T a t is future AT at time t month; T w t is the predicted water temperature value at time t month.
For the analyses, we focused on the catchment's observed data periods (Fig. 1b) and future periods 2021-2050 and 2071-2100, followed by the 30 years for a climatological standard normal 53 .
Oxygen saturation. Waters with concentrations below saturation are called "deficit" whereas those with concentrations exceeding saturation are called "supersaturated". As a result, the oxygen saturation concentration serves as the baseline for any endeavor to measure oxygen-based water quality by determining the oxygen concentration of unpolluted water 1 . The saturated DO concentration depends on the temperature, salinity of water, and oxygen partial pressure. Saturated DO concentration is influenced by these elements, as indicated by 54 where o s = saturated DO concentration (mgO 2 /L), ω k , ω s = elevation above sea level (dimensionless), and salinity (dimensionless) respectively, and o sf = the saturated DO concentration of sea-level freshwater (mgO 2 /L). The following are the individual impacts of temperature, salinity, and elevation. where T abs = absolute temperature in kelvin.
Salinity, S (ppt). The oxygen saturation of seawater is calculated by multiplying the sea-level freshwater saturation by 54 Elevation, k (km). The influence of atmospheric pressure on gas saturation at elevation is based on the standard atmosphere as described by the cubic polynomial 54 Additional insight into DO can be obtained by computing the rate of change of saturation by differentiating Eq. (2) with respect to temperature. Although functions like Eq. (2) can sometimes be differentiated analytically, the results are cumbersome and typically provide no insight. Numerical differentiation provides an alternative means to obtain the same results with the centered divided difference 55 where x = the value of the independent variable, h ′ (x) = the function's first derivative with respect to x evaluated at x, and = a very small perturbation of x. For the present case, with x = T and h(x) = o s (T) , the result is do s (T)/dT with units of (mgO 2 /L)/°C.    Fig. 1a, and their main characteristics with study periods are outlined in Fig. 1b in tabular form. Two data sources were used to compile the models, with one being global, and one regional. We have used the Global Freshwater Quality Database (GEMSTAT) data for Narmada, Cauvery, Sabarmati, and Godavari catchments. We have used the Central Water Commission (CWC), India data for Tunga-Bhadra, Musi, and Ganga catchments. The Global Freshwater Quality Database GEMStat 56 is hosted by the International Centre for Water Resources and Global Change (ICWRGC) and provides inland water quality data within the framework of the GEMS/Water Programme of the United Nations Environment Programme (UNEP). Approximately 500 water quality parameters were available in the global GEMSTAT database, out of which water temperature was used in this study for Narmada, Cauvery, Sabarmati, Godavari catchments when compiling models. The gauging stations are run by the Central Water Commission (CWC), India, and measure water temperature (T w ) over a period of time (monthly mean of ten samples) 57 . We observed that majority of the time series retrieved from the source datasets (GEMSTAT and CWC) are discontinuous. To build a kNN-LSTM model, a complete dataset is necessary. To build an entire data record, the na.interp() method in R's forecast library was utilized to interpolate the missing observations using the STL (Seasonal and Trend decomposition using Loess) decomposition 58 . The meteorological data used in this work are monthly minimum (T min ), and maximum (T max ) air temperatures. T min , T max was available from the India Meteorological Department (IMD) data on a 1° Latitude × 1° Longitude grids spatial resolution from 1951 to 2018. We have spatially interpolated the AT observations to the RWT gauging locations using linear interpolation. We averaged T min and T max to get the monthly mean AT as widely used literature 59 . Figure 1b shows the catchment means for all variables. This study used the subset of the National Aeronautics Space Administration (NASA) Earth Exchange Global Daily Downscaled Projections (NEX-GDDP) dataset to assess the impact of climate change on RWTs for seven catchments of India. The NEX-GDDP is made up of downscaled climate scenarios for the entire world produced from the General Circulation Model (GCM) runs undertaken as part of the Coupled Model Intercomparison Project Phase 5 (CMIP5) and spanning two of the four greenhouse gas emissions scenarios known as RCPs 60 . The ensemble mean of the NEX-GDDP dataset contains RCP 4.5 and RCP 8.5 downscaled projections from the 21 GCMs models and scenarios, and each climate projection has daily maximum temperature, minimum temperature, and precipitation for 1950 through 2100. The dataset has a spatial resolution of 0.25° (~ 25 km × 25 km). This study retrieved the daily T min and T max values, converted them into a monthly scale, and averaged them to obtain the monthly mean AT for future RWT predictions.
Data pre-processing. The applied data pre-processing consists of aggregating multiple data sources and feature engineering. We examined the data's autocorrelation and partial autocorrelation functions (ACF and PACF) to account for the time-lag information in RWT prediction at monthly time scale. These functions suggest that the 1-month time-lag is significant in the observed record. Thus, air temperature (AT[t]), and time-lag effects of air and water temperatures (AT[t − 1], RWT[t − 1]) are used as input variables in the prediction of RWT.
In the k-NN bootstrap resampling algorithm, "one simulation" is defined as a set of simulated values with a length equal to the observed dataset and chosen to generate 50 simulations. Following that, we ran a comparison study of monthly statistics (maximum, minimum, mean, standard deviation) for both the historical and simulated ensemble records. Also, we compared the lag-1 autocorrelation of the k-NN simulated data with observed data. The comparison has revealed that the algorithm produced the applicable distributional statistics of the observed dataset, implying that the algorithm generates accurate and diverse conditions. The lag-1 autocorrelation represents the relationship between two consecutive time steps (e.g., x t and x t-1 ). When we compare the lag-1 autocorrelation of the historical and simulated record, we find that the lag-1 autocorrelation's seasonality is frequently reproduced. Then, the simulated values of the monthly average ATs and RWTs were then used as model input in the LSTM model. In this study, while training a kNN-LSTM model on a time series, all the possible combinations of LSTM hyperparameter sets (the number of LSTM hidden layers: 1-3, the total number of units per layer: 5-100, time steps:1-12, the dropout ratio: 0-0.4, epochs: 50-100, and the batch size: 2-64) are evaluated using an emerging state-of-the-art Bayesian Optimization approach to optimize the hyperparameters, and the topmost group is chosen to improve the model's performance.

Results
The data used in this work comprises monthly average AT and the corresponding RWT for seven majorly polluted river locations in India. We used meteorological definitions of seasons: monsoon = June, July, August, September; post-monsoon = October, November; winter = December, January, February; and summer = March, April, May 61 . The catchment means of RWT, and AT for all seven catchments ranged between 24.68 °C, 30.34 °C, and 24.24 °C, 28.81 °C, respectively (Fig. 1b).
To examine the variability of annually averaged AT, RWT, and DO changes, the study calculated the linear trends using the observed data for seven catchments of India (Fig. 2). The AT and RWT increased and observed DO has decreased during the studied period for all catchments except Cauvery, Godavari, and Ganga catchments (Fig. 2). The RWT rising rates are lower than those of AT in general.
Air temperature has shown a rising trend except for Cauvery (− 0.01 °C/year) catchment, and the rising rates range from 0.002 to 0.380 °C/year. RWT shows a rising trend except for Cauvery (− 0.06 °C/year), Godavari (− 0.03 °C/year), and Ganga (− 0.07 °C/year) catchments, and the rising rates vary between 0.01 and 0. 17  DO shows a decreasing trend except for Cauvery (0.01 (mgO 2 /L)/year), Godavari (0.004 (mgO 2 /L)/year), and Ganga (0.01 (mgO 2 /L)/year) catchments (where there is a significant decreasing trend of AT and RWT has been noted), and the decreasing rates vary between − 0.01 and − 0.003 (mgO 2 /L)/year. Such DO decrease patterns have been explored in several locations throughout the world. The DO, for instance, has been a seasonal DO variation, low (DO < 10 mgO 2 /L) and high (DO > 14 mgO 2 /L) over Clackamas River near Oregon City, OR, USA 14 , and rising RWTs in the Delaware River, the USA by 2 °C to peak summer levels of 30 °C, based on saturation, DO levels will decline by about 0.2 mgO 2 /L 13 . Generally, RWT and AT are directly correlated, but RWT and DO are inversely correlated 14,16 . However, for the Godavari and Ganga catchment, the water temperature has shown decreasing trend (− 0.03 °C/year and − 0.07 °C/year respectively) with an increasing trend of AT (0.01 °C/year and 0.08 °C/year, respectively), which specifies that the temporal shifts of RWT may not be explained AT alone. RWT is directly influenced by multiple parameters, including streamflow 68 , river geometry, groundwater inputs, slope, water depth, etc. 69 .
Deep learning model performance. The simulated samples from the k-NN bootstrap resampling algorithm and lag variables as input to the kNN-LSTM hybrid model to predict the RWT for monthly data for all the seven catchments of India. To mathematically quantify the predictive performances of kNN-LSTM approach, six statistical measures are calculated, such as the coefficient of determination (R 2 ), Kling-Gupta efficiency (KGE) 70 , RMSE-observations standard deviation ratio (RSR) 71 , the root mean squared error (RMSE), Nash-Sutcliffe efficiency (NSE) 72 , and the mean absolute error (MAE) (Fig. 3). Detailed descriptions of these metrics can be found in Rajesh et al. 35 .
The relationship between monthly RWT and AT at seven catchments is relatively strongly correlated for the kNN-LSTM model (R 2 and NSE values). The RMSE metrics varied from 1.266 to 2.361 for kNN-LSTM monthly data estimated between observed and simulated for all the catchments (Fig. 3). The NSE values for all the catchments range from 0.446 to 0.920 (Fig. 3) for the kNN-LSTM model for monthly data, which is reasonable compared with earlier standalone LSTM models by Stajkowski et al. 73 73 used AT values as input for hourly data in their analysis, Qiu et al. 74 used AT and discharge as input for daily data in RWT predictions, and the current study is dedicated to monthly timescales. Based on RSR, KGE, R 2 , and NSE performance values (Fig. 3), the kNN-LSTM model is the best performant model for all catchments. Overall, the kNN-LSTM model statistical metrics are reasonably within the range for all the catchment locations providing confidence that the developed model performs effectively.
The following analyses concentrate on how RWT affects the oxygen saturation of seven catchments of India.   Table 1). The magnitude of DO decrease with respect to minimum RWT increase is higher for Narmada, Sabarmati, Godavari, and Ganga catchments, and with respect to maximum RWT increase is higher for Narmada, Musi, and Ganga catchments (Table 1). Overall, results indicated that DO with respect to RWT over Indian catchments would likely drop by more than 0.02 (mgO 2 /L)/°C for 2071-2100 (Table 1). Figure 5a shows the rate of change of DO saturation levels under mean river water temperature do s (T)/dT ((mgO 2 /L)/°C) for observed and projected (2071-2100) data for seven Indian catchments. The vertical dotted lines indicate the mean of historical (Tw hist °C) and projected (2071-2100) (Tw proj °C) water temperatures. As depicted in Fig. 5a, projected (2071-2100) (Tw proj °C) water temperatures increase is higher for Narmada, Tunga-Bhadra, Musi, and Ganga catchments compared to historical (Tw hist °C), which leads to a higher drop in the rate of change of oxygen saturation for these catchments. For Cauvery catchment, projected (2071-2100) (Tw proj °C) water temperatures increase is low compared to historical (Tw hist °C), which leads to a minimal drop in the rate of change of oxygen saturation. www.nature.com/scientificreports/ The specification of a DO water-quality standard, o wq (mgO 2 /L), is used to evaluate oxygen assimilative capacity. Figure 5b shows the DO concentration (mgO 2 /L) scale with respect to the observed (blue color) and projected (2071-2100) (red color) minimum, maximum and average water temperature (°C) levels of seven Indian catchments. From Fig. 5b, observed that 10.3, 6.6, and 7.9 mgO 2 /L, and 9.1, 6.3, and 7.3 mgO 2 /L DO concentrations with respect to historical and projected (2071-2100) minimum, maximum and average water temperatures (°C) respectively of seven Indian catchments. DO concentration (mgO 2 /L) scale scores are dropped from 7.9 to 7.3 mgO 2 /L respective to the observed and projected (2071-2100) mean RWT levels of seven catchments (Fig. 5b). Table 2 listed the DO concentrations and DO decrease percentage with respect to monthly average summer and winter RWTs for historical and projected (2071-2100) with RCP 8.5 experiments for seven Indian catchments. The summer RWT increase for Tunga-Bhadra, Sabarmati, Musi, and Ganga basins are predicted as 3.1, 3.8, 5.8, 7.3 °C, respectively, with a more pronounced increase of 7.8 °C for the Narmada River for 2071-2100. The magnitude of DO concentrations decreases with respect to summer RWT increases is higher for Narmada, Musi, and Ganga catchment sites, and the percentage of DO decreases for 2071-2100 relative to the historical values noted 12.4, 9.3, and 11.9%, respectively ( Table 2). The low DO concentrations decrease was observed for Cauvery and Godavari, and the percentage of DO decrease was noted as 1.0 and 3.3%, respectively (Table 2).  www.nature.com/scientificreports/ Overall, the summer displayed larger percent decreases in DO compared to the winter season, and the largest DO decreases were found in the Narmada catchment (Table 2).

Discussion
This study presents new intuitions on the assessment of climate change impacts on saturated DO concentrations with respect to RWT for seven different catchment sites across India in different physiographic settings. For this, using the monthly kNN-LSTM prediction model, which is developed based on AT, including time-lag effects, demonstrates rising RWTs will reduce a river's assimilative capacity by affecting its oxygen metabolism, in addition to lowering saturation. The monthly R 2 scores estimated between observed and simulated RWT using kNN-LSTM for various stations ranged between 0.446 and 0.920, KGE scores ranged between 0.378 and 0.868, NSE scores ranged between 0.446 and 0.920, RSR scores ranged between 0.283 and 0.744, and RMSE scores were ≤ 2.4 °C during the testing periods for kNN-LSTM prediction model, revealing high model reliability. All the developed model statistical metrics covered the range of model reliability described in the literature. The RMSE scores for all the catchments ranged between 1.266 and 2.361 °C pertaining to the kNN-LSTM model for monthly data, which are reasonable in comparison to earlier models of the Spatio-temporal approach by Jackson et al. 75 (1.570 °C); Bayesian regression approach by Sohrabi et al. 68 73,74,76 . However, it can be noted that the study was conducted by Feigl et al. 76 used AT, runoff, precipitation, and global radiation values as input in the RWT prediction for daily data, the study by Qiu et al. 74 used daily AT, and discharge as input in RWT prediction, and the study by Stajkowski et al. 73 used AT values as input in RWT prediction for hourly data.
The RWT increases of up to7 °C for summer, reaching close to 35 °C, decreases DO by 2-12%, thus decreasing the saturation capacity for DO for 2071-2100. DO concentration (mgO 2 /L) scale scores are dropped from 7.9 to 7.3 mgO 2 /L respective to the observed and projected (2071-2100) mean RWT levels of seven catchments. These scores reveal that DO concentration (mgO 2 /L) values are dropping for projected years as RWTs rise. The RWT increases of up to7 °C for summer, demonstrated in this work, were found to agree with Chapra et al. 1 (5 °C increments in summer RWTs in most of the world's rivers over the next 50 years). The DO concentration (mgO 2 /L) scale scores, as demonstrated in this work, were found to agree with Du et al. 15 (DO concentrations on the basin average scale will decrease by 0.72 mgO 2 /L under RCP 8.5 scenario for 2061-2100) and Chapra et al. 1 (DO oxygen concentrations are 9.0 and 6.8 mgO 2 /L for freshwater temperatures 20 and 35 °C, respectively).
The percentage of DO decrease with respect to summer RWTs is higher for Narmada, Musi, and Ganga catchment sites for 2071-2100 relative to the historical values noted as 12.4, 9.3, and 11.9%, respectively, probably because of the influence of disposal of untreated sewage and industrial wastewater along with due to increased reaction kinetics at a higher temperature under climate change scenarios (Table 2). In this study, overall, for all seven catchments, the decrease in DO is 8% for the plausible future (2071-2100) ( Fig. 5b and Table 2). These projected change patterns are most consistent with earlier hydrological model studies by Ficklin et al. 16 (10% decreases in DO by 2100 at Sierra Nevada in California, USA) and by Du et al. 15 (DO decrease on the basin average scale by 0.72 mgO 2 /L under RCP 8.5 scenario for 2061-2100 in the Athabasca River Basin, Canada). Overall, this study demonstrated how river oxygen levels would be influenced by rising RWT due to climate change using the kNN-LSTM model for the Indian riverine system. The rising RWTs will reduce river assimilative capacity by affecting its oxygen metabolism, in addition to lowering saturation, and necessitates redefining/alterations of the river water quality standards under climate change.
Furthermore, the DO simulated by Eq. (2) is the saturated oxygen concentration, which is the total amount of DO that can be dissolved within the streamflow volume, and thus, it can be expected that the DO concentrations presented in this study represent the ceiling of potential DO levels. Though the hybrid kNN-LSTM model performed well, further research is needed to improve it. We found that inherent uncertainties from the kNN-LSTM model can accumulate and affect the final performance measurements. Such uncertainties can originate from various sources, from noise and temporal discontinuity present in the original water quality sampled observations to the model hyperparameters used to predict RWT. To address such model uncertainties, the ensemble of DL models, etc. 79,80 can be adopted which can combine RWT predictions from various DL models, allowing the decision-makers to choose the best possible prediction within a range of predictions 81 . Despite the effectiveness of the modeling frameworks, as demonstrated in the present work, it has some limitations. Firstly, flow discharge may play a vital role in RWT predictions, especially the Indian rivers are significantly impacted by low flows in the summer season. However, flow discharge was not examined in this work due to a lack of complete data. Therefore, we will strengthen the hybrid modeling framework in future research by integrating flow discharge as model input for rivers. Secondly, RWT is directly influenced by multiple parameters, including streamflow 28,65,68 , river geometry, groundwater inputs, slope, water depth, etc. 69  www.nature.com/scientificreports/ in the present study. This study set the Salinity (S) value for seven river catchments to zero because most rivers and streams had minimal salinity 1 . Overall, this research offers vital intuitions about the historical and projected RWT and DO states of major Indian river catchment locations, which may be beneficial in creating future water management plans that may impact aquatic resources.

Conclusions
The study demonstrates the climate change impacts on saturated DO concentrations with respect to RWT for the seven major polluted Indian catchments at a monthly timescale. The hybrid kNN-LSTM model is implemented in this study to predict the RWT addressing sparse spatiotemporal RWT data. Further assessed the climate change impacts on DO concentrations with respect to RWT using a forced by an ensemble of RCP 8.5 scenario downscaled projections of AT data from the NEX-GDDP dataset. The results lead to the following conclusions: 1. An increase in AT will have an effect on RWTs, and saturated DO concentrations. The latter will trigger higher RWT and lower DO concentration. These changes appear especially significant for the summer seasons and include RWT increases of up to7 °C for summer, reaching close to 35 °C, decreases of DO by 2-12%, thus decreasing the saturation capacity for DO. 2. The percentage of decrease of DO saturation levels with respect to summer RWTs is higher for Narmada, Musi, and Ganga catchment sites for 2071-2100 relative to the historical values noted as 12.4, 9.3, and 11.9%, respectively. 3. DO concentration (mgO 2 /L) scale scores are dropped from 7.9 to 7.3 mgO 2 /L respective to the observed and projected (2071-2100) mean RWT levels of seven catchments. 4. Overall, saturated DO concentration (mgO 2 /L) levels are dropping by 8% under the rise of summer RWT by more than 4.3 °C for 2071-2100. That is, for every 1 °C RWT increase, there will be about 2.3% decrease in DO saturation level concentrations over Indian catchments under climate signals. 5. The study provides an assessment of the individual contribution of RWT rise on depletion of saturated DO levels, which is helpful for the policymakers and pollution control authorities for sustainable river water quality management. www.nature.com/scientificreports/